Realization of higher Wess-Zumino-Witten models in spin chains 
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Building on the generalization of the exactly dimerized Majumdar-Ghosh ground state to arbitrary 
spin S for the Heisenberg chain with a three-site term (Sj_i ■ Sj)(Sj ■ Si+i) + H.c, we use density- 
matrix renormalization group simulations and exact diagonalizations to determine the nature of the 
dimerization transition for S = 1, 3/2 and 2. The resulting central charge and critical exponent are 
in good agreement with the SU(2)t=2S Wess-Zumino-Witten values c = 3fc/(2 + k) and r\ = 3/(2 + 
k). Since the 3-site term that induces dimerization appears naturally if exchange interactions are 
calculated beyond second order, these results suggest that SU(2)k>i Wess-Zumino-Witten models 
might finally be realized in actual spin chains. 



INTRODUCTION 

The solution in 1931 of the one-dimensional spin- 
1/2 Heisenberg model by Bethe[l] is the first example 
of a long series of exact solutions of integrable models 
that includes models of direct physical relevance such as 
the one-dimensional Hubbard model[2] and the Kondo 
problem[3, 4]. As far as spin chains are concerned, there 
exist two families of integrable spin-S models which can 
be viewed as generalizations of Bethe's solution. In both 
cases, the coupling is limited to nearest neighbours, and 
the coupling between spins at sites i and i + 1 takes the 
form of a polynomial of degree 2S in Si ■ S i+1 . In the first 
family, the couplings are such that the Hamiltonian is the 
sum of permutations between nearest neighbours, so that 
the symmetry is enlarged to SU(N) with N = 25 + 1[5]. 
For this family, the central charge is equal to 25*, and 
the low-energy theory is equivalent to 2S free bosons. 
Another family of integrable models has been discov- 
ered in the early eighties in which the ground state is 
also critical[6, 7], but the low energy field theory is the 
SU(2)k=2S Wess-Zumino-Witten model [8], with central 
charge c = 3fc/(2 + k)[9]. These two classes of model 
have been later on shown to belong to a single, more gen- 
eral family of integrable SU(N) models with higher sym- 
metric representations[10, 11], with a low-energy sector 
described by the SU(N)k Wess-Zumino-Witten model, 
where k is the rank of the symmetric irreducible repre- 
sentation of SU(N), as shown with Bethe ansatz[12] and 
confirmed numerically [13]. In both classes, the Hamilto- 
nian contains higher powers of Sj ■ S^+i with significant 
coefficients, and the possibility to realize them in actual 
spin chains is quite remote. For the first family of mod- 
els, an alternative to spin chains is provided by ultracold 
fermionic alcaline earth atoms [14], for which the SU(N) 
Heisenberg model is the relevant effective model in the 
Mott insulating phase, but for the second family, there is 
no such alternative. 

Shortly after the discovery of the second family of 
integrable models described by the SU(2)i,=2S Wess- 



Zumino-Witten universality class, it has been suggested 
that the simple spin-3/2 Heisenberg chain, which is ex- 
pected to be critical, might fall into this category[15], but 
this proposal has been definitely excluded by numerical 
simulations[16], which have shown that the model has 
a central charge equal to 1. Since then, other models 
whose low-energy physics is described by the SU(2) k=2S 
Wess-Zumino-Witten universality class have been discov- 
ered, including a generalization of the Haldane-Shastry 
model to arbitrary spins[17], but all these models con- 
tain spin-spin interactions that seem impossible to realize 
in actual systems, and the search for spin chain models 
that are both realistic and described at low energy by 
the SU(2)k=2S Wess-Zumino-Witten universality class is 
still open. 

Another route to this universality class, not followed 
by any concrete implementation so far, has been out- 
lined by Affleck and Haldane[8], who argued on the basis 
of renormalization group arguments that the transition 
from a uniform to a dimerized phase might take place 
through a critical point of universality class 5 l J7(2)j, =1 
or 5f/(2)fc =2 5 depending on the initial conditions of the 
flow. If the model is explicitly dimerized, the transition 
is expected to be in the usual SU (2)k=i class, but if there 
is no explicit dimerization, as in the bilinear-biquadratic 
spin-1 chain, then the critical point might be SU (2) k=2S ■ 

In that respect, the recent generalization to arbitrary 
spin S of the Majumdar-Ghosh exactly dimerized ground 
state of the J\ — J2 spin-1/2 chain, where J\ and J 2 are 
the exchange constants between nearest and next-nearest 
neighbors[18, 19], opens new perspectives. Indeed, it has 
been shown that the spin-S" Heisenberg Hamiltonian with 
a three-site term 

n = ^(JiS 4 -s 4+1 + j 3 [(Si • s i+ i) (s i+1 • s J+2 ) + h.c.}) 

i 

(i) 

is equivalent, for S — 1/2, to the Ji — J2 model with 
J 2 = J3/2, and that, for arbitrary S, the ground state is 
exactly dimerized when J 3 / Ji = 1/(45(5 + 1) — 2). Since 
the ground state of the pure Heisenberg model ( J3 = G) is 
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not dimerized, there has to be a dimerization transition 
upon increasing J3, and, since the model is not explic- 
itly dimerized, the transition may have the more exotic 
SU(2)k = 2s class for arbitrary 5, as already shown for 
5=1 [20]. 

In this paper, using extensive density matrix renor- 
malization group (DMRG) [21-23] simulations and exact 
diagonalizations (ED), we show that this is indeed the 
case for 5 = 3/2 and 5 = 2. First of all, using a combi- 
nation of level spectroscopy [24-27] and DMRG results 
for the size dependence of dimerization, we precisely lo- 
cate the transition to dimerization, with the conclusion 
that it occurs for quite small ratios of the coupling J3/J1 
(0.063 and 0.0403 respectively). Then, using DMRG es- 
timates of the central charge (extracted from the entan- 
glement entropy) and of the spin-spin correlations, we 
show that, at the critical point, and to a very good ac- 
curacy, the central charge is equal to c = 3fc/(2 + k) 
and the exponent of the spin-spin correlation function is 
equal to 77 = 3/(2 + k), with k = 25, in agreement with 
the SU(2)k=2S WZW universality class. So the model of 
Eq. (1) appears to be a realistic one to implement the 
SU{2)k=2S WZW universality class in actual spin chains. 

Numerical methods. — The main results of this paper 
have been obtained using DMRG simulations. To com- 
pute the correlation functions and the dimerization, wc 
have used open boundary conditions. We have performed 
12 sweeps keeping up to 1600 states in the last sweep. 
This has allowed us to reach good convergence up to 150 
sites, with a discarded weight of the order of 10~ 8 close 
to the transition (and smaller far from the transition). 

The calculation of the central charge is more difficult 
because it requires calculations with periodic boundary 
conditions, where DMRG is less efficient. We have de- 
duced the central charge from calculations of the entan- 
glement entropy on systems with L = 30, 40, 50 and 60 
sites for 5 = 3/2, and L = 30, 40 and 50 sites for 5 = 2. 
We did at least 7 sweeps, keeping 3500 states in the fi- 
nal sweep. Even with this high number of states, the 
maximum discarded weight was about 10 -6 . To get ac- 
curate estimates, we had to perform an extrapolation of 
the central charge with respect to the discarded weight, 
as explained below. 

Exact diagonalizations of finite clusters have also been 
used in the context of the so-called level spectroscopy 
method[24-27] to locate the transition. To that end, we 
have used the Lanczos algorithm to compute the ground 
state energy in different sectors of total magnetization 
Sf 0t , momentum k and parity under the transformation 
i — > L — i, where % is the site label and L the total number 
of sites. For spin 5 = 3/2 chain, we have diagonalized 
systems with up to L = 16, and for spin 5 = 2 chains, 
we went up to L = 12. 

Determination of the critical point. — We have used 
two different techniques to get an accurate estimate of 
the location of the transition between the uniform and 



the dimerized phase, one based on a direct calculation of 
the order parameter associated to the dimerization with 
DMRG, the other one based on exact diagonalizations. 

To determine the location of the phase transition with 
DMRG, we have computed the order parameter associ- 
ated to the dimerization d = \Si-\ ■ Si — Si ■ Si+i\. The 
order parameter in the center of the chain for different 
system sizes is shown in Fig. 1 for spin 5 = 3/2 system, 
and in Fig. 2 for spin 5 = 2. First of all, the data on 
finite-size systems are in both cases consistent with a con- 
tinuous transition. As for the spin-1/2 Jj — J 2 chain, the 
fact that the spin-3/2 chain is gapless below the transi- 
tion, hence that the spin stiffness jumps at the transition, 
does not affect the continuous nature of the dimerization 
transition. In the thermodynamic limit, the dimerization 
order parameter is thus expected to be zero in the uni- 
form phase, and to start growing at the transition with a 
critical exponent /3 : d cx ((J3 — Js c )/Js)^ ■ However, it 
is not possible to locate the transition precisely using this 
scaling property because of finite-size effects. It is much 
more efficient to rely on the fact that, at the transition, 
the order parameter should tend to zero as a power law 
as a function of the system size, while it should tend to a 
finite value above it. Therefore, to obtain the transition 
point, we have plotted d as a function of the system size L 
on a log-log scale. At the transition, we expect a straight 
line, while above the transition we expect a convex be- 
haviour. This analysis is done in the left inset of Fig. 1 
for spin 5 = 3/2 and of Fig. 2 for spin 5 = 2. It leads 
to a rather precise determination of the critical points at 
J Zc jJ x = 0.063 for 5 = 3/2 and at J 3c /Ji = 0.0403 for 
5 = 2. Note that below the transition, the behaviour is 
expected to depend on the value of the spin: for inte- 
ger spins, the spectrum is gapped, and the dimerization 
should decay exponentially, while for half-integer spin, 
the system is gapless with a Luttinger liquid exponent 
K = 1, and the dimerization should decay as the inverse 
of the square of the size. Interestingly, the behaviours 
below the transition look very similar for spin 3/2 and 2, 
showing the same kind of concavity. For spin 3/2, this 
is probably due to the fact that the fast decay only ap- 
pears for larger sizes, so that for the sizes available one 
only observes a cross-over regime, while for spin 2, the 
concavity reflects the exponential decay. 

To check the reliability of these results, we have per- 
formed, following Ref. 24-27, a level spectroscopy analy- 
sis. Since the transition is second order, we do not expect 
any level crossing in the ground state. However, the na- 
ture of the first excited state is expected to change at the 
transition [24] . For a system with L — 4n, in the uniform 
region, the first excited state is a triplet with wave vec- 
tor 7r and parity P — — 1, while in the dimerized region, 
the first excited state is a singlet with again wave vector 
7r but P = 1. This method has been shown to lead to 
very accurate results [24] even if only rather small sizes 
are available because the finite-size scaling of the crossing 
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Figure 1. (Color online) Dimerization as a function of J3/J1 
for S = 3/2 for different system sizes with up to L = 150 sites 
in the vicinity of the phase transition J3/J1 ~ 0.063. The left 
insets show the size dependence at J3/J1 = 0.0625, 0.063 and 
0.0635, respectively. The right inset shows the localisation 
of the first excited state crossing point as a function of the 
system size. 
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Figure 2. (Color online) Dimerization as a function of J3/J1 
for 5* = 2 for different system sizes with up to L = 150 sites in 
the vicinity of the phase transition J3/J1 « 0.0403. The left 
insets show the size dependence at J3/J1 = 0.0402, 0.0403 
and 0.0404, respectively. The right inset shows the localisa- 
tion of the first excited state crossing point as a function of 
the system size. 



point is very fast: a c (L) = a c + constant/L 2 + 0(1/L 3 ), 
i.e. deviations scale as the square of the inverse size with- 
out logarithmic corrections. 

For spin 5 = 3/2, we have calculated the first excited 
states for three system sizes (L = 8, 12 and 16). The 
results are plotted in the right inset of Fig. 1. The scal- 
ing with 1/L 2 is quite accurate, and the level crossing 
extrapolates to J 3c / ' J\ = 0.064, in good agreement with 
the DMRG result J 3c /Ji = 0.063. For spin S = 2, we 
have calculated the first excited states for two system 
sizes (L = 8, L = 12). The results are shown in the 
right inset of Fig. 2. The level crossing extrapolates 



to J 3c /Ji = 0.0406. The 1/L 2 behavior could not be 
checked in this case because we have only two points. The 
result is nonetheless in good agreement with the DMRG 
result J 3c / J 1 = 0.0403. 

Central charge. — At the transition between the uni- 
form and the dimerized phase, it has already been shown 
that the 5 = 1 chain with Heisenberg and three-body 
interactions can be effectively described by a SU(2) k= 2S 
WZW model[20]. To check if this result can be extended 
to larger 5, we have calculated numerically two quanti- 
ties which are fixed by the SU(2) k= 2S WZW universal- 
ity class, namely the central charge, which is given by 
c = 3k/ (2 + k), and the critical exponent of the decay 
of the spin-spin correlation function, which is given by 
7] = 3/(2 + *). 

The central charge can be extracted for a finite system 
from the entanglement entropy obtained with DMRG on 
a system with PBC. More precisely, it is obtained by 
fitting the numerical results with the following formula 
[281: 



5, = -In 



~ L ■ 1 




— sin 
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where c is the central charge, Si — — Ti'QelnQe the en- 
tanglement entropy, and ge the reduced density matrix 
of a subsystem of size I. This analysis should then be 
followed by a finite size scaling of c. However, due to the 
rapid increase of the cpu time required, we have not been 
able to reach full convergence at the critical point for 30 
sites and beyond. For a given size, we have therefore kept 
track of the discarded weight at various steps of the sim- 
ulation, and we have performed an extrapolation of the 
central charge with respect to the discarded weight. A 
linear function turned out to fit well the results as shown 
in the inset of Fig 3 for S = 3/2 and in the inset of Fig. 

4 for 5 = 2, leading to estimates of the central charge 
for different sizes. This analysis has been performed for 
L = 30,40,50,60 for 5 = 3/2 and L = 30,40,50 for 

5 = 2, and a finite-size scaling of the result has revealed 
a linear behaviour as a function of 1 /L. For spin 5 = 3/2, 
the extrapolation to the thermodynamic limit is shown 
in Fig. 3. The SU(2) k=3 WZW theory predicts a cen- 
tral charge of 9/5 = 1.8. Our numerical extrapolation 
leads to c = 1.807, which is remarkably close to the ex- 
pected value, leaving little doubt about the nature of the 
transition. For 5 = 2, the extrapolation shown in Fig. 
4 leads to c = 2.01, again in very good agreement with 
the value c = 2 for the SU(2) k=i WZW model. These 
results nicely confirm that the transition is also in the 
SU{2) k=2S WZW universality class for S=3/2 and S=2. 

Critical exponent of the spin-spin correlation 
function. — We now turn to the extraction of the 
critical exponent associated to the spin-spin correlation 
function defined by C s (i,j) — (5f ■ 5|). At the transi- 
tion, it is expected to decay as a power law according to 
C s (i,j) oc (i — j) v . In such a case, the boundary effects 



4 



1.88 




0.01 0.02 0.03 0.04 0.05 



l/L 

Figure 3. (Color online) Extrapolation to the thermodynamic 
limit of the central charge for S = 3/2. The data points 
shown are results from an extrapolation as a function of the 
discarded weight, which is displayed in the inset. 
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Figure 5. (Color online) Decay of the spin correlation function 
(5f ■ Sj) with distance \i — j\ for systems with L = 150 and 
spin 5 = 3/2 and 5 = 2. The solid lines are fits to power laws 
with exponents indicated in the plot. 
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Figure 4. (Color online) The same as in Fig. 3 for S=2. 

are expected to play an important role, and one can 
only hope to observe the bulk behavior not too close to 
the boundaries. In practice, we could only perform a fit 
to determine rj inside a narrow window in the centre of 
the chain where the decay looked polynomial. In the 
fit, we only considered odd distances to get rid of the 
inaccuracy coming from the odd-even oscillation in the 
correlation function. The extrapolations are shown in 
Fig. 5. For 5 = 3/2, we have obtained rj = 0.62. This 
has to be compared to the WZW model where rj = 3/5. 
So, within a small region, the correlation function seems 
to be well described by the WZW model. The same 
analysis for 5 = 2 leads to an exponent rj — 0.52, which 
should be compared to the result 77 = 1/2 of the WZW 
theory. These results are therefore compatible with a 
SU{2) k=2S WZW theory. 

Discussion and conclusion. — To summarize, we have 
shown that, for 5 = 3/2 and 5 = 2, the dimeriza- 
tion transition in the J\ — J3 model between a uni- 
form and a spontaneously dimerized phase is in the 



SU{2)k = 2s WZW universality class. The evidence re- 
lies on a careful extrapolation of the central charge at 
the transition, which agrees with the theoretical predic- 
tion c = 3k/ (2 + k) with an error smaller than 1 per- 
cent (see Table I) . To further confirm this conclusion, we 
have computed the exponent r\ of the spin-spin correla- 
tion function, and the results are again compatible with 
a WZW critical point. All the numerical results obtained 
in this paper, and in a former paper [20] for 5 = 1, are 
summarized in Table I. 

In view of these results, it is tempting to conjecture 
that the universality class of the dimerization transition 
will be SU(2)k=2S for all values of 5. This is likely to 
be the case if the transition is continuous. However, 
as pointed out in the context of one-dimensional cou- 
pled dipolar gases, the fixed point might be unstable for 
5 > 2 because of the presence of relevant operators that 
are forbidden by symmetry for 5 < 2 [29], leading to a 
first order transition, in agreement with recent numeri- 
cal results[30]. This issue, which requires to investigate 
spins larger than 2, is left for future investigation. 





5=1 


5 = 3/2 


5 = 2 


25 = k 


J 3c /Ji (DMRG) 


0.111 


0.063 


0.0403 




J 3c /Ji (ED) 




0.064 


0.0406 




J3/J1 (MG point) 


1/6 


1/13 


1/22 


1/(45(5 4- 1) - 2) 


c (WZW) 


3/2 


9/5 


2 


3k/(2 + k) 


c (DMRG) 


1.502 


1.807 


2.01 




n (WZW) 


3/4 


3/5 


1/2 


3/(2 + k) 


■q (DMRG) 


0.72 


0.62 


0.52 





Table I. Summary of all quantities computed with DMRG 
and the same quantities in the WZW model: the ratio J3/J1 
at the dimerization transition and at the MG point, the cen- 
tral charge c and the exponent rj of the spin-spin correlation 
function. The data for S=l are from Ref. [20]. 
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Finally, by contrast to the spin S = 1/2 J\ — J2 chain 
where the critical point is the standard SU{2) level 1 
field theory, the physics is not standard for higher spin. 
This is particularly clear for spin S = 3/2, where the 
central charge takes the rather exotic value c = 9/5. To 
our knowledge, up to now no simple and realistic spin 
chain model which allows to probe this kind of physics 
has been put forward. In that respect, the Hamiltonian 
we propose has many advantages: it is simple - it has 
only two different terms and interactions do not reach 
beyond next-nearest neighbors - and it is realistic since 
the three-site term emerges naturally from a strong cou- 
pling expansion of a multi-orbital Hubbard model. It is 
thus our hope that the results of the present paper will 
prove useful on the way to an experimental realization of 
higher WZW models in spin chains. 

We acknowledge useful discussions with D. Cabra, F. 
Essler and T. Ziman. We are especially indebted to P. 
Lcchcminant for several interesting comments, in partic- 
ular for pointing out to us the possibility of a first order 
transition for spins larger than 2. This work has been 
supported by the Swiss National Fund and by MaNEP. 
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